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(N ! ABSTRACT 

Pure luminosity evolution models for galaxies provide an unacceptable fit to the redshifts and 
colors of faint galaxies. In this paper we demonstrate, using HST morphological number counts 
derived both from the /si4"band of WFPC2 in the Medium Deep Survey (MDS) and the Hubble 
Deep Field (HDF) and from the i2i.6-band of NICMOS, and ground-based spectroscopic data 
of the Hawaii Deep Field and the Canada-France Redshift Survey, that number evolution is 
necessary for galaxies, regardless of whether the cosmic geometry is flat, open, or A-dominated. 
Furthermore, we show that the number evolution is small at redshifts of z < 1, but large at 
| z > 1, and that this conclusion is valid for all the three cosmological models under consideration. 

If the universe is open or A-dominated, the models, which are subject to the constraint of the 
conservation of the comoving mass density of galaxies, naturally predict a population of star- 
forming galaxies with the redshift distribution peaking at z = 2 ~ 3, which seems to be 
consistent with the recent findings from Lyman-break photometric selection techniques. If the 
cosmological model is flat, however, the conservation of the comoving mass density is invalid. 
Hence, in order to account for the steep slope of -B-band number counts at faint magnitudes in 
the flat universe, such a star-forming galaxy population has to be introduced ad hoc into the 
> ■ modelling alongside the merger assumption. 
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Telescope (HST), whose unprecedented imaging abil- 
ity enables galaxies to be segregated morphologically 
into several broad classes (Glazebrook et al. 1995a; 
Driver et al. 1995; Abraham et al. 1996; Teplitz et 
al. 1998). These morphological data have provided 
us with valuable information to understand the origin 
and evolution of galaxies. Now, any self-consistent 
theory of galaxy evolution and cosmology must pass 
the test of matching the distribution of galaxy prop- 
erties derived from HST morphological data. 

In a previous investigation (He & Zhang 1998, here- 
after HZ98), we modelled the number counts of E/SO 
galaxies obtained from the MDS and the HDF in the 
HST /gu-bandpass (A e jj=8000A), and found that 
the number counts of ellipticals could be well ex- 
plained by PLE models in any cosmological model 
if ellipticals are assumed to have formed at high 
redshifts (say Zf=5.0) and thereafter have passively 
evolved (i.e., no further star formation). This is just 
the traditional scenario of the formation and evolu- 
tion for elliptical galaxies. However, incorporating 
other observations such as redshift and color distribu- 
tions and the modelling of the other types (i.e., early- 
and late- type spirals and irregulars, as separated by 
HST), we have found that the FBG problem cannot 
be resolved by the conventional scenario for ellipti- 
cal galaxies, and the inclusion of dust extinction also 
cannot save the scenario. This conclusion holds in 
each of the three cosmological models under consider- 
ation: flat, open and A-dominated. Our investigation 
suggests that number evolution may be essential to 
explain the faint data (He & Zhang 1999a, hereafter 
HZ99a). 

Number evolution models have been widely consid- 
ered to account for the paradoxes concerning the FBG 
problem, whereas PLE model is the starting point 
of our investigation. We model the number counts 
of the three classes of galaxies- ellipticals, early- type 
spirals, and late- type spirals/irregulars- which are de- 
rived from the HST Igu-passband, and in combina- 
tion with the analysis of redshift distributions in the 
observed magnitude bins, we find that PLE models 
can explain the number counts at the bright magni- 
tudes of Isu < 22.5. This indicates that PLE models 
could be valid up to z ~ 1, in apparent contradic- 
tion to scenarios of strong number evolution at low 
redshifts (e.g., Kauffmann, Chariot, & White 1996). 
Kauffmann et al. concluded that only one third the 
number of nearby bright elliptical galaxies were in 
place at z ~ 1, which would seem to be inconsistent 
with the number counts at the bright end. Our con- 
clusion is consistent with those of Roukema & Yoshii 
(1993) and Woods & Fahlman (1997) based on the 
analysis of angular correlation function, and that of 
Totani & Yoshii (1998) based on the V/V max anal- 



ysis (Schmidt 1968), and in particular, agrees well 
with the results of Im et al. (1999), based on the 
morphologically divided redshift distribution of faint 
galaxies, and those of Brinchmann et al. (1998) and 
Lilly et al. (1998), based on the HST imaging of the 
CFRS and LDSS redshift surveys. 

We tentatively construct a set of phenomenologi- 
cal number-luminosity evolution (NLE) models to ex- 
plain data such as the morphological number counts 
in both the Isi4-band of the Wide Field and Plane- 
tary Camera 2 (WFPC2) and the Fi. 6 -band of the 
Near-Infrared Camera and Multiobject Spectrometer 
(NICMOS), as well as the counts in B- and X-bands 
derived from most of the ground-based observations 
and redshift distributions derived from the Hawaii 
sample (Cowie et al. 1996) and the Canada-France 
Redshift Survey (CFRS) (Lilly et al. 1995a). The 
significant feature of these models is that the num- 
ber evolution is small at redshifts of z < 1, while it 
is large at high redshifts. We find that NLE models 
can roughly reproduce these observations. 

To simplify the modelling, we employ the con- 
straint of conservation of the comoving mass den- 
sity of galaxies, in the form of Guiderdoni & Rocca- 
Volmerange (1991) . We find that a population of late- 
type spiral/irregular galaxies and star- forming galax- 
ies naturally emerge in our models, due to mass con- 
servation. Such a population will peak at z ss 2.8 — 3. 
Between z ~ 3 and the present day, it behaves as 
a population which seems to be fading, disappearing 
or self-destructing, with very blue colours. Such a 
population also seems to be consistent with the re- 
cent findings using Lyman-break selection techniques 
(see Steidel & Hamilton 1992, 1993; Steidel, Pettini, 
&: Hamilton 1995; Madau 1995). This conclusion, 
however, is only valid in open (e.g., rio=0.1) or A- 
dominated (e.g., f2o=0.2, Ao=0.8) cosmological mod- 
els. If the universe is flat (Oo=l), the comoving mass 
density of such a predicted population becomes nega- 
tive at high redshifts (z > 3), which is physically un- 
reasonable. Therefore, to account for the faint data, 
we abandon the constraint of mass conservation and 
introduce a star-forming population ad hoc into the 
modelling. 

In Section 2, we present the PLE models. We are 
particularly interested in the behaviors of the mod- 
els at bright magnitudes, or at low redshifts. Section 
3 is devoted to the number evolution models. The 
summary and conclusions are presented in Section 4. 

2. PURE LUMINOSITY EVOLUTION 

PLE model is the starting point of the investiga- 
tion. This kind of model takes as ingredients the 
stellar evolutionary tracks, normally for a restricted 
metallicity range, the initial mass function (IMF), 
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and a gas consumption timescale associated with star 
formation rate (SFR), which are adjusted to give the 
present range of colors across the Hubble sequence. 
With each galaxy assumed to evolve as an isolated 
system, the rest-frame spectral energy distribution 
(SED) can be predicted at a given time and thus 
an evolutionary correction can be determined (Ellis 
1997). The evolutionary SEDs of galaxies are com- 
puted using population synthesis models, (e.g., the 
latest Bruzual & Chariot synthesis code, which we 
employ in this work). The reader is referred to HZ98, 
HZ99a or Pozzetti et al. (1996) for details of the mod- 
elling. We also consider the effects of dust extinction, 
following Wang's (1991) prescription, which is simi- 
lar to that of HZ98 and HZ99a. However we choose 
a larger optical depth of r* = 0.3 than that adopted 
by Wang for galaxies other than ellipticals (see also 
Campos & Shanks 1997). 

Similar to the prescription of HZ98 and HZ99a, 
we adopt three representative cosmological models in 
this work: 1) flat, 1\) = 1.0, Ao = 0, and h = 0.5 
(#0=100 h km s _1 Mpc -1 ) (hereafter, Scenario A), 
2) open, fio = 0.1, Ao = 0, and h = 0.5 (Scenario 
B), and 3) flat and A-dominated, £\) = 0.2, Ao = 0.8, 
and h = 0.6 (Scenario C). We assume a formation 
redshift Zf = 5.0 for all types of galaxies in all three 
cosmological models. 

We first consider the modelling of the morpholog- 
ical number counts derived from the MDS and the 
HDF in the /g^-passband. The basic ingredient of 
the models, the present-day luminosity function (LF), 
is not well-determined by local surveys due to in- 
herent uncertainties such as a large local fluctuation 
(Shanks 1989), significant local evolution (Maddox 
et al. 1990a), selection effects and/or incompleteness 
(Zwicky 1957; Disney 1976; Ferguson & McGaugh 
1995), or perhaps systematic errors in local surveys 
(Metcalfe, Fong & Shanks 1995a). Therefore, a sim- 
ple 3 parameter Schechter function ((/>*, M* and a, 
cf., Schechter 1976), as adopted by most authors, 
seems not to be universal, and hence not to be quite 
appropriate as an average over all types and environ- 
ments (cf., Binggeli, Sandage, & Tammann 1988). In 
particular, much recent work finds a tail for dwarfs 
and late types that is steeper than a ~-1.3 (Driver, 
Windhorst & Griffiths 1995), and in addition, low 
surface brightness galaxies (LSBGs) may make up a 
significant amount of the luminosity density of the lo- 
cal universe (see Impey &: Bothun 1997 for a review 
about LSBGs). Regardless of these uncertainties, the 
simple Schechter form for local LFs is still our work- 
ing assumption in the current research. 

Firstly, we adjust the Schechter LF parameters by 
trial and error so that the model predictions can 
match better the bright-end number counts, since the 



count normalization is sensitive to the three param- 
eters (f)*, M* and a, especially the former two (cf. 
Ellis 1997). The LFs for galaxies are listed in Table 
1. Secondly, we test these ad hoc adoptions by com- 
parison with observations. It can be seen, in Figure 

1, that both the B- and the K- band adopted LFs 
lie close to the observations, and in particular, the 
observed faint-end slopes are not as high as, say -1.5, 
as predicted by the semi-analytic models (cf. Cole et 
al. 1994; Kauffmann, White & Guiderdoni 1993). 

2.1. Number counts for elliptical galaxies 

We adopt the Scab (1986) IMF for ellipticals and 
the Salpeter (1955) IMF for the other types. As ad- 
dressed in HZ98, the Scalo IMF is less rich in massive 
stars than the Salpeter IMF because of the steeper 
slope of the former at the high-mass end. The SFR 
for ellipticals is adopted as the simple exponential 
form 

Y>(t) ~ exp(-t/r e ), (1) 

where ip(t) refers to the SFR, and r e is the timescale 
characterizing this kind of SFR. A single value of r e 
can not satisfactorily reproduce the B — K color dis- 
tribution (He & Zhang 1999b), which may be a re- 
flection of intrinsic complexities such as the metal- 
licity variations among the population of ellipticals. 
Besides, there is also an age-metallicity degeneracy 
(Worthey 1994), which greatly complicates the sit- 
uation. For the sake of simplicity, we neglect the 
chemical evolution in the current work. In our mod- 
els, ellipticals are subdivided into five parts equal in 
number, with each one of them specified by a differ- 
ent r e for SFR to characterize its luminosity evolu- 
tion. By such treatments, we hope to circumvent the 
metallicity-age degeneracy. The values of r e , mod- 
elled B — I and B — K colors are shown in Table 

2. Notice that the observed average B-K color for 
present-day ellipticals is ~ 4.10 (Johnson photomet- 
ric system), with the dispersion being about 0.1 — 0.2 
mag, and hence these values of the parameter r e are 
acceptable. It is well known that the local photo- 
metric and spectroscopic properties of ellipticals are 
largely insensitive to the past star formation history 
(cf. Bruzual 1983; Bower, Lucey & Ellis 1992; Char- 
lot & Silk 1994). 

We show in Fig.2-(a) the predicted number counts 
for ellipticals by PLE models to compare with those 
derived from HST in ig^-passband. It can be seen 
that up to /8i4=22.5, the predictions can satisfacto- 
rily reproduce the bright-end number counts in all 
the three scenarios, indicating the PLE models are 
valid within this magnitude range. At the faint-end 
(^814 > 22.5), however, the predictions in Scenarios 
B and C largely overestimate the observations. Al- 
though the prediction for Scenario A fits the bright 
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counts well, as we have analyzed in HZ99a, models 
with strong luminosity evolution predict too many 
high-z objects, in substantial disagreement with the 
color-selected sample from Cowie et al. (1996). 

Though PLE models are unacceptable, we present 
redshift distributions for ellipticals in the observed 
I8i4 magnitude bins, in order to better understand 
the models. In Fig.2-(b), (c) and (d), we show pre- 
dictions of redshift distributions in the three mag- 
nitude bins, I 8 u < 21.0, 21.0 < l 8 u < 22.5, and 
22.5 < ^814 < 25.0. The redshift distributions peak 
around z ~ 0.5, indicating the counts at Isu < 21.0 
are contributed by the objects at the redshifts of 
z < 0.5 or even as large as z ~ 0.8 (the half maximum 
of the peaks), and from Fig.2-(c), it can be seen that 
the peaks appear at z ~ 0.8, indicating the objects 
between 21.0 < Isu < 22.5 are mostly at z < 0.8 or 
~ 1.0. From the previous analysis, we know that the 
number counts at Isu < 22.5 can be well explained 
by PLE models up to z ~ 1. As addressed above, the 
predicted number counts at Isu > 22.5 in Scenarios 
B and C significantly exceed the observations, and 
from Fig.2-(d), we know that these excessive objects 
mostly lie in z > 1. This implies, for example, that if 
a merger scenario for the evolution of ellipticals is in- 
volved, then the merger rate in the interval < z < 1 
seems to be small. This conclusion is valid in all the 
three world models under consideration. 

We now consider an example of strong number evo- 
lution for ellipticals. Kauffmann, Chariot & White 
(1996) claimed, based on < V/V max > analysis 
(Schmidt 1968) and the data from the CFRS and 
the Hawaii Deep Survey (Cowie et al. 1996), that 
only one third of the nearby bright elliptical galaxies 
existed at z ~ 1. We adopt their expression of the 
number density evolution for ellipticals with respect 
to z, i.e., F = (1 + z)~ 7 , with 7 = 1.5 ± 0.4 (This 
means Fpsl/3atz = l), to compare the predictions 
with observations. We show the model predictions 
in Fig.2-(a). Obviously, such a strong number evolu- 
tion cannot be accepted, since it underestimates the 
bright-end number counts. Notice that the bright- 
end counts are mostly contributed by L* galaxies, 
and hence they are insensitive to the LF faint-end 
slope a. An enhancement of the normalization <p* by, 
say 50% can improve the goodness of fit, but such 
an enhancement will produce too many bright ellip- 
ticals at the present day, and conflict with the local 
surveys, especially the LF at If -band. It is worthy 
of mentioning that Totani & Yoshii (1998) have per- 
formed a detailed V/V max test, and concluded that 
the value of < V/V max > for ellipticals is ~0.5 in the 
redshift range 0.3 < z < 0.8, which means that the 
number density is nearly unchanged in this redshift 
range, and this conclusion is similar to ours. We con- 



clude that strong number evolution for ellipticals at 
low redshifts is not acceptable. 

2.2. Number counts for early-type spiral galaxies 

Using the morphological classification of galaxies 
by HST, we combine Sab and Sbc galaxies into a 
group of early-type spirals. The LFs are presented 
in Table 1. The SFR is taken the same exponen- 
tial form as Eq. [1], but the timescales are generally 
larger than those of ellipticals in order to match their 
local photometric properties. These parameters are 
listed in Table 2. 

In Fig.3-(a), we present our PLE model predictions 
of number counts to compare with the observations. 
The models satisfactorily reproduce the observations 
at almost all magnitudes in Scenarios B and C, indi- 
cating little number evolution. For Scenario A, how- 
ever, the PLE model cannot account for observations 
at igu > 22.5, and hence other assumptions than 
PLE are needed. As is Section 2.1, we analyze the 
redshift distributions in the observed magnitude bins. 

We can see from Fig.3-(b) and (c) that, in the mag- 
nitude bins Isu < 21.0 and 21.0 < I 8 u < 22.5, there 
are peaks in the redshift distributions at z ~ 0.5 ~ 

0. 8, indicating that PLE models could be valid even 
out to the redshift of the half maximum of the peaks, 

1. e., z ps 0.8 ~ 1.0. For I 8 u > 22.5, we can see 
from Fig.3-(d) that the models fail to predict enough 
early-type spirals at the redshifts of z > 1.0. 

2.3. Number counts for late-type spiral/irregular 

galaxies 

We mix Scd and Sdm/Irr classes into a category of 
late- type spiral/irregular galaxies. In addition, we 
introduce another population which has very blue 
colors (even bluer than irregular galaxies), and we 
refer to it as the vB galaxies with the name taken 
after Pozzetti et al. (1996). The vB galaxies are in- 
cluded to explain the bluest colours, B — K ps 2 of 
the Hawaii sample and (V — I)ab ~ 0.5 of the CFRS 
sample. As addressed by Totani & Yoshii (1998), 
such very blue galaxies comprise nearly 10% of the 
CFRS sample. The very blue colors suggest that 
they are star-forming galaxies, and hence vB galaxies 
are meant to reproduce the population of star-burst 
galaxies present at each z, whose evolution does not 
follow the prescription of standard population syn- 
thesis. Following the treatment of Pozzetti et al., we 
assume that star formation in these galaxies keeps 
their SEDs constant in time, and the SEDs can be 
approximated by a model with constant SFR lasting 
1.5 Gyr. 

Their LFs are also listed in Table 1. In Fig.4-(a), 
we show the predicted number counts for comparison 
with observations. Similar to the cases of ellipticals 
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and early-type spirals, the PLE models can also re- 
produce the late-type spiral/irregular population to 
1814 < 22.5, indicating that PLE models could be 
valid in the redshift range z < 0.8 ~ 1 (Fig.4-(b) and 
(c)), while at Isu > 22.5, especially for Scenario A, 
the models underpredict the observations, indicating 
the models need to be improved at high redshifts of 
z > 1 (Fig.4-(d)). 

In short, PLE models can well explain the number 
counts of all types of galaxies at the magnitudes of 
Isu < 22.5, or equivalently, according to the analy- 
sis of redshift distributions, at the redshifts of z < 1. 
Hence, strong number evolution at low redshifts is 
not acceptable, comparing to the observations. 

3. NUMBER-LUMINOSITY EVOLUTION 

Pure luminosity evolution works fairly well at low 
redshifts, z < 1, as suggested by our previous in- 
vestigations. Next, without considering the underly- 
ing physical mechanism, we phenomenologically con- 
struct a set of number evolution models to explain 
the faint data. There is no doubt that the luminos- 
ity evolution must occur, due to the aging of star 
populations. Moreover, if mergers between galax- 
ies are included in the modelling, it is more real- 
istic to consider the additional star formation, that 
can be induced by interactions between galaxies. For 
simplicity, we use the same galaxy evolution mod- 
els as in Section 2 to compute the X-corrections and 
evolutionary-corrections for galaxies, without con- 
sidering more complicated star formation scenarios 
(e.g., Colin, Schramm, & Peimbert 1994; Fritze- 
v.Alvensleben & Gerhard 1994), or various compli- 
cated physical processes concerning the formation of 
galaxies (cf. Roukema et al. 1997). 

Number evolution can be simply realized by mak- 
ing the LF parameters </>*, M* and a functions of z. 
For this, we adopt the following expressions: 

<f (z) = <fo(l + ^) Ql,+Q < (2) 

where 4>q refers to the characteristic density of 
present-day galaxies. Qo and Q2 are two free parame- 
ters. This form is a refinement of our previous formal- 
ism (cf. HZ99a), and can also be treated as a mod- 
ification of that of Guiderdoni & Rocca-Volmerange 
(1991), in which the power-law index is not a func- 
tion of z. Obviously, Q2 will mostly affect the high-z 
behavior. The evolution of the faint-end slope a is 
parametrized as the following: 

a(z) =00-71 {£) ~ 72 (^) ■ (3) 

where cto represents the local quantity, and zj is the 
formation redshift. The two parameters 71 and 72 



characterize the functional relationship between a 
and z. In principle, M* should also be a function 
of z. We find, however, the variation of M* can be to 
some extent compensated by that of (f>* and a, and 
for simplicity, we neglect the variation of the charac- 
teristic magnitude. Our model is therefore similar to 
the 4>*-a model of Guiderdoni & Rocca-Volmerange 
(1991). 

Thus we have four free parameters, Qo, Q2, 71 and 
72, with which to characterize the number evolution. 
Nevertheless, for a specific Hubble type, there is no 
need to invoke simultaneously all the four parame- 
ters. Our strategy to proceed is to involve as few of 
these parameters as possible, as long as the models 
can reproduce the observed data. We model the three 
Hubble types independently. 

3.1. Morphological number counts 

We first consider the explanation of elliptical num- 
ber counts by the NLE models. From Section 2, we 
realize that the number evolution is small at z < 1, 
hence the (absolute) value of the parameter Qo should 
be low. Since the PLE predictions exceed the obser- 
vations at faint magnitudes (i.e, at high redshifts), 
the parameter Q2 should be negative. We try to fix 
these two parameters by trial and error to match the 
data. It can be seen from Panel (a) of Figure 5, that 
the NLE models can reproduce well the data at both 
bright and faint magnitudes in each of the three Sce- 
narios. The values of Qo and Q2 are listed in Table 3, 
and in Figure 6, we show the relationship of 0* with 
z. We can see that the number density of ellipticals is 
nearly unchanged in the redshift range of z < 1; while 
at high redshifts, ellipticals are absent regardless of 
the cosmological models. This conclusion agrees well 
with the finding of Totani & Yoshii (1998) based on 
the V/V max analysis, and is also in agreement with 
that of Zepf (1997). 

The number counts of early-type spirals can also be 
satisfactorily reproduced by our models (See Panel 
(b) of Figure 5). For Scenarios B and C, only 71 
and 72 are needed, and Qo or Q2 is not necessary. 
There are additional advantages for such a-variation 
models, that the faint-end counts are contributed by 
dwarfs rather than bright galaxies, and hence there 
are no high- z tails predicted by the models. For Sce- 
nario A, it seems that Q2 is also needed, since the flat 
universe has less volume at high redshifts than the 
open or A-dominated universe, and hence a slightly 
positive Q2 can help to improve the fit at faint magni- 
tudes. The values of these parameters are also listed 
in Table 3. We can see that the variation of a at low 
redshifts is also not dramatically large. For exam- 
ple, a steepens only by 0.06 from z = to z = 1 in 
Scenario A. 
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The modelling for the late- type spiral/irregular and 
vB galaxies is similar to that for ellipticals and early- 
type spirals, but requires even more free parameters. 
As a result, we proceed in an alternative way. In a 
scenario of hierarchical clustering, galaxies are con- 
stantly merging over cosmic time. If the hierarchical 
picture is realistic, and in the absence of other scenar- 
ios such as fading or self-destructing of a galaxy pop- 
ulation, then the total comoving masses of gas and 
stars can be regarded to be enclosed within galax- 
ies. Hence the total comoving mass density should 
be conserved in merger processes. We employ the 
conservation of mass to constrain the number den- 
sity evolution of the latest Hubble types. 

Mass (baryonic, including gas and stars) conserva- 
tion can be expressed in the following way: 

M*(z)(f)*(z)r(a j (z) + 2) = constant, (4) 

3 

which is taken from Guiderdoni & Rocca-Volmerange 
(1991). M* is the characteristic mass of the galaxy 
mass function (MF). The total mass density is com- 
puted by summing over galaxy type j, and the 
constant can be determined by present-day quanti- 
ties. In order to obtain the evolution of LFs, we 
need to convert the evolving MFs into evolving LFs, 
and hence the mass-to-luminosity (M/L) ratios are 
needed. Unfortunately, the M/L ratios are poorly 
determined, especially for £?-passband. However, the 
near-infrared light is less affected by star formation, 
and thus by evolution, than the optical emission, and 
it is a good tracer of the total mass within galaxies. 
Hence we will make use of if-band LFs instead of the 
MFs of galaxies. 

We assume the faint-end slope a for late-type 
galaxies (Scd, Sdm/Irr, and vB) evolves in the same 
way as for the early- type spirals (Sab and Sbc). In 
this case the evolution of (j)* can be determined by 
Eq. [4]. The evolution of the characteristic number 
density (j>* for L* galaxies is illustrated in Figure 7. 
We first analyze the cases for Scenarios B and C. Us- 
ing the constraint of mass conservation suggested by 
Eq. [4], (f>* increases with increasing z to the max- 
imum at z ~ 2.8 in both Scenarios B and C, and 
then decreases to almost zero at the formation red- 
shift Zf. At z = 1, the number density (/)* for L* 
galaxies is 44% and 88% higher than the present-day 
value for Scenarios B and C, respectively. Notice that 
the Sdm/Irr/vB galaxies are star- forming galaxies, 
as suggested by their blue colors and predicted by 
our population synthesis models. If we only focus on 
these galaxies, they appear to be a disappearing or 
self-destructing population (cf. Babul & Rees 1992; 
Babul & Ferguson 1996; Ferguson h Babul 1998), but 
they are produced naturally by the merger scenario, 



and are not needed to be ad hoc introduced into the 
models. In particular, since the faint-end slope a also 
increases with z, there are even more dwarf galax- 
ies at high redshifts. Such a population seems to be 
qualitatively consistent with the recent findings us- 
ing Lyman-break selection techniques (see Steidel & 
Hamilton 1992, 1993; Steidel, Pettini, & Hamilton 
1995; Madau 1995) that a substantial population of 
star-forming galaxies exists at z > 3 (Steidel et al. 
1996; Steidel et al. 1998). We believe that such a 
population of Lyman-break galaxies can be identified 
with the Sdm/Irr /vB galaxies suggested by our mod- 
els. 

Such an analysis, however, cannot be applied to the 
case of Scenario A. In Figure 7, <j>* becomes negative 
when z > 3, which is obviously unphysical. The rea- 
son why <p becomes negative at high redshifts is that 
a flat universe has less volume at a given redshift than 
an open or A-dominated universe. Therefore, in or- 
der to account for the faint-end counts, the number 
evolution of early-type spirals for Scenario A has to 
be larger than that for Scenarios B and C (See Table 
3). On the other hand, the number evolution of ellip- 
ticals for Scenario A is less than for Scenarios B and 
C. As a result, the increase of the mass density for 
early-type spirals exceeds the corresponding decrease 
for ellipticals. Therefore, if the model is subject to 
the constraint of the mass density conservation, it 
will inevitably predict a negative <p* for Sdm/Irr/vB 
galaxies at high-redshifts. The implication is that 
the mass density conservation suggested by Eq. [4] 
cannot be applied to the case of Scenario A, and fur- 
thermore, in order to explain the faint data, the mass 
density must increase with z. 

The failure of the assumption of mass density con- 
servation in Scenario A can be understood in several 
different ways, and its cosmological interpretation is 
also not unique. There may be a disappearing/self- 
destructing population of star-forming galaxies at 
high redshift which has no present-day counterpart. 
Another possibility is that galaxies have been losing 
mass during their evolution, in particular, during the 
period of merging between galaxies, due to super- 
novae which enable the gas to escape from the poten- 
tial well of dark matter halos. The situations may be 
very complicated. 

Regardless of the difficulty in explaining the phys- 
ical mechanism behind the non-conservation of the 
mass density for galaxies, we can follow the same ad 
hoc prescription as that for early-type spirals to ad- 
just the parameters (Qo, Q2, 71, and 72) of the models 
to match the data, without discriminating between 
the two possibilities discussed above. The values of 
these parameters are also listed in Table 3. The pre- 
dicted number counts for the latest type galaxies are 
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shown in Panel (c) of Figure 5. It can be seen that the 
predictions can be accepted within the uncertainties. 
In addition, the predictions for the overall population 
(shown in Panel (d) of Figure 5) are also acceptable. 

The morphological number counts in the near- 
infrared bands have also been provided by NICMOS 
(Teplitz et al. 1998). Predictions of the NLE mod- 
els are presented in Figure 8 to compare with the 
observed data in the Hi q band (A e //= 1.6/xm). We 
can see that the models can reasonably reproduce the 
data. 

3.2. Number counts in B- and K- passbands 

We now examine whether the number evolution 
models can reproduce B- and K-bands number 
counts obtained from the majority of surveys so far. 
From Figure 9 we can see that the models can re- 
produce better the number counts at bright magni- 
tudes in both the B (B < 21) and the K (K < 17) 
bands, indicating that the normalization (involving 
both (J)q and L* , cf. Ellis 1997) is appropriate using 
the adopted LFs. (Due to large uncertainties in the 
local LFs, the normalizations are usually chosen to 
scale the model predictions to the observed number 
counts at a relatively faint magnitude of B ~ 19.0. 
See Pozzetti et al. 1996 for a comprehensive re- 
view). Furthermore, the NLE models can reproduce 
the counts at the faintest B and K magnitudes for 
each of the three Scenarios. 

3.3. Redshift distribution 

The spectroscopic studies with the LRIS spectro- 
graph on the Keck Telescope of two of the Hawaii 
deep survey fields SSA 13 and SSA 22 (Cowie et al. 
1996) and the CFRS present two of the largest and 
deepest redshift samples so far, providing powerful 
tools for understanding the evolution of galaxies. 

In Fig.lO-(a), (b), and (c), we show the predicted z- 
distributions in the magnitude range of 22.5 < B < 
24.0 in the three Scenarios, to compare with mea- 
sured redshifts from Cowie et al. (1996), whose com- 
pleteness is more than 85%. We can see that the 
predictions roughly match the observed data, and in 
particular, there are no very high-z tails. In Fig. 10- 
(d), (e), and (f), we compare the NLE-model pre- 
dicted z-distributions, limited in the magnitude range 
of 17.5 < I A b < 22.5, with the CFRS sample (Lilly 
et al. 1995a; Le Fevre et al. 1995; Lilly et al. 1995b; 
Hammer et al. 1995; Crampton et al. 1995). Notice 
that the CFRS data are less deep than the Hawaii 
sample, with the completeness ~85%. We can see 
that the predictions are reasonably matched by the 
data, especially for Scenarios B and C. However, we 
should also realize that all models overpredict galax- 
ies with z > 1. The overprediction may be due to 



the z-unidentified galaxies, since the completeness of 
both of the two samples is around 85%, and, accord- 
ing to Pozzetti et al. , such z-unidentified galaxies are 
most probably at high z (see Section 3.6 of Pozzetti 
et al. 1996). 

4. SUMMARY AND CONCLUSIONS 

Recent morphological data on faint galaxies are of 
great significance for us in understanding the dis- 
tant universe, and have shed new light on the long- 
standing FBG problem. We have systematically in- 
vestigated this question using HST data. In previous 
work, we found that the number counts of ellipticals 
could be explained successfully by PLE models only 
a scenario where ellipticals have extremely short SFR 
timescales (r e = 0.2, 0, and 0.05 Gyr for flat, open 
and A -dominated universes, respectively), i.e., they 
formed at high-z (say zj = 5) and have passively 
evolved since their formation (HZ98). Such a view, 
however, would contradict other observations such as 
color distributions (HZ99a). Hence, the evolution of 
galaxy number density seems to be required, and one 
of the main purposes of this work is to study the na- 
ture of the number evolution. 

PLE models are the starting point of our investi- 
gation. By modelling the number counts of the three 
morphological classes in the HST /si4-passband, in 
combination with the analysis of redshift distribu- 
tions in the three observed magnitude bins, isi4 < 
21.0, 21.0 < hiA < 22.5, and 22.5 < I 8 u < 25.0, we 
find that PLE models can fit the number counts at 
the bright magnitudes of Isu < 22.5, indicating such 
PLE models could be valid up to z « 1. As a result, 
scenarios of strong number evolution at low redshifts 
such as that of Kauffmann et al. (1996) would seem 
to be rejected. The results agree well with those of Im 
et al. (1999), based on the morphologically divided 
redshift distribution of faint galaxies, and those of 
Brinchmann et al. (1998) and Lilly et al. (1998), 
based on the HST imaging of the CFRS and LDSS 
redshift surveys. In addition, our conclusions are also 
consistent with those of Roukema & Yoshii (1993) 
and Woods & Fahlman (1997) based on analysis of 
the angular correlation function, and that of Totani 
& Yoshii (1998) based on the V/V max analysis. 

We then constructed a set of phenomenological 
number-luminosity evolution models and by trial and 
error fixed a set of model parameters to explain the 
data such as morphological number counts both in 
I8i4- and in .Big-bands, as well as the counts in the 
B- and i^-bands, and redshift distributions derived 
from the Hawaii and the CFRS samples. The sig- 
nificant feature of these models is that the number 
evolution is small at redshifts of z < 1 and large at 
higher redshifts. We find that the models can roughly 



reproduce the above-mentioned observations in each 
of the cosmologies under consideration. 

To simplify the modelling, we employ the con- 
straint of conservation of the comoving mass den- 
sity of galaxies, in the form presented by Guiderdoni 
& Rocca-Volmerange (1991). We find that a popu- 
lation of late- type spiral/irregular galaxies and very 
blue galaxies naturally emerge in these models. Such 
a population will peak at z ~ 2.8. We argue that 
such a population seems to be consistent with the 
recent findings using Lyman-break photometric se- 
lection techniques. This conclusion, however, is only 
valid in an open or A-dominated cosmological mod- 
els. If the universe is flat, the comoving number den- 
sity of such a predicted population becomes negative 
at high redshifts [z > 3), and hence it is obviously 
physically unacceptable. Therefore, in order to ac- 
count for the faint data, such a population has to be 
introduced ad hoc into the modelling. The failure of 
the mass conservation assumption in a flat universe 
model, as discussed in Section 3.1, may imply the 
existence of a disappearing/self-destructing popula- 
tion of star- forming galaxies at high redshifts, whose 
present-day counterpart is not detectable. Alterna- 
tively, galaxies may have been losing mass during 
their lifetimes, due to supernova bursts, for exam- 



ple. If so, the evolutionary scenario for galaxies may 
be extremely complicated. 

We have not taken into account the evolution of 
metallicity with redshift, and details such as the recy- 
cling of the residual gas ejected by dying stars. These 
would be future improvements over our current sim- 
ple prescriptions. Nevertheless, the models faithfully 
reproduce the number counts and the redshift distri- 
butions, and these results can be treated as heuristic 
guidelines for further investigations. 
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Table 1 

-b-band lf parameters for present-day galaxies. 





E/SO 


Sab 


Sbc 


Scd 


Sdm/Irr 


vB 




10.0 


4.00 


6.00 


6.00 


4.00 


4.00 


M* 


-21.0 


-20.80 


-20.80 


-20.80 


-20.70 


-20.50 


a 


-0.90 


-0.95 


-1.10 


-1.10 


-1.20 


-1.30 



Notes: 

(1) . <fi* , M*, and a are determined by trial and error, see text for details; 

(2) . (f>* in units of 10~ 4 /Mpc 3 ; 

(3) . vB galaxies, taken after Pozzetti et al. (1996), are included to explain 
the bluest samples, e.g., CFRS. See text for details; 

(4) . The Hubble constant has been scaled to h=0.5. 



Table 2 

SFR PARAMETER T e (MEASURED IN Gyr) AND MODELLED B — I AND B-K COLORS 
FOR GALAXIES IN THE THREE COSMOLOGICAL MODELS UNDER CONSIDERATION. 



Type 




Seen. 


A 




Seen. 


B 




Seen. C 






B - I 


B-K 


Te 


B - I 


B-K 


T e 


B-I 


B-K 


E/S0-I 


0.5 


2.29 


4.14 


1.0 


2.34 


4.17 


1.1 


2.35 


4.18 


E/S0-II 


0.7 


2.29 


4.12 


1.4 


2.33 


4.16 


1.6 


2.33 


4.16 


E/S0-III 


1.0 


2.28 


4.11 


2.1 


2.30 


4.11 


2.4 


2.30 


4.10 


E/S0-IV 


1.5 


2.25 


4.07 


3.0 


2.24 


4.06 


3.1 


2.24 


4.07 


E/S0-V 


2.1 


2.21 


4.00 


4.1 


2.15 


3.94 


4.3 


2.15 


3.94 


Sab 


2.3 


2.15 


3.89 


3.3 


2.17 


3.93 


3.5 


2.15 


3.91 


Sbc 


4.6 


1.81 


3.45 


6.3 


1.84 


3.51 


7.0 


1.81 


3.47 


Scd 


9.5 


1.57 


3.14 


13.5 


1.60 


3.19 


15.0 


1.59 


3.18 


Sdm/Irr 


oo 


1.37 


2.87 


oo 


1.42 


2.95 


oo 


1.43 


2.96 


vB 




0.92 


2.25 




0.92 


2.25 




0.92 


2.25 



Notes: 

(1) . We choose a Scalo IMF for ellipticals, and a Salpeter IMF for the other types; 

(2) . T e = oo for Sdm/Irr galaxies denotes constant SFR; 

(3) . The SFR of vB galaxies is not described by Eq. [1]. See text for detail. 



Table 3 

The phenomenological parameters for the number evolution models in the 
three cosmological models under consideration. 



Parameter 


E/SO 


Sab 


Sbc 


Scd 


Sdm/Irr 


vB 


Seen. A 


Qo 


-0.02 











1.0 


1.0 


Q 2 


-0.05 


0.03 


0.03 


0.03 


0.03 


0.03 


7i 





0.2 


0.2 


0.2 


0.2 


0.2 


72 





0.5 


0.5 


0.5 


0.5 


0.5 


Seen. B 


Qo 


-0.05 















Q 2 


-0.15 















7i 





0.2 


0.2 


0.2 


0.2 


0.2 


72 





0.4 


0.4 


0.4 


0.4 


0.4 


Seen. C 


Qo 


-0.2 















Q 2 


-0.15 















71 





0.2 


0.2 


0.2 


0.2 


0.2 


72 





0.4 


0.4 


0.4 


0.4 


0.4 



Notes: 

Qo and Q2 for Sdm/Irr and vB galaxies in Scenarios B and C are determined by Eq. [4]. 
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Fig. 1. — Present-day luminosity functions of galaxies. The adopted LFs are indicated by thick lines, whose values correspond 
to those in Table 1. Here and hereafter, the capital letters A, B and C refer to Scenarios A, B and C, respectively. Observational 
data were obtained as follows, (a) Solid line with dots: Efstathiou, Ellis, & Peterson (1988), dotted line with triangles: Loveday 
ct al. (1992), dashed line with squares: Marzke et al. (1994), dashed line with stars: Zucca ct al. (1997). (b) Solid line with 
dots: Gardner et al. (1997), dotted line with triangles: Glazebrook et al. (1995b), dashed line with squares: Szokoly et al. (1998; 
C-method), dashed line with stars: Szokoly et al. (1998; conventional method). Notice that the characteristic density <f>* and the 
mixing ratio between different galaxies of Efstathiou et al. (1988) LF are taken from Pozzetti et al. (1996). Panels (a) and (b) are 
for B and K bands, respectively. The Hubble constant has been scaled to h = 0.5 for Scenario C. 
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Fig. 2. — (a). Differential number counts for E/SO galaxies as a function of apparent magnitude in Isi4 band. The sources of the 
observational data are exhibited in the figure and these data are indicated by symbols. The model predictions are shown by lines. 
'PLE' denotes pure luminosity evolution models, and 'SNE', predictions of the strong number evolution, based on the formalism 
of Kauffmann et al. (1996). See Section 2.1 for details. Fig.2-(b), (c), and (d) are the z distributions for E/SO limited to three 
magnitude bins, I$i4 < 21.0, 21.0 < 7~8i4 < 22.5 and 22.5 < /8i4 < 25.0, respectively. Models are indicated by lines, and the model 
predictions have been normalized in units of deg~ 2 with a redshift bin Az = 0.01. 
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Fig. 3. — The same as Figure 2, but for early-type spiral galaxies (Sab+Sbc). 
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Fig. 4.— 



The same as Figure 2 and 3, but for late-type spiral/irregular (Scd+Sdm/Irr) and vB galaxies. 
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(a) Ellipticals 

■k MDS , Glazebrook et al. (1995a) 
A MDS, Driver etal. (1995) 
• HDF, Abraham etal. (1996) 



(b) Early-type Spirals 

* MDS , Glazebrook et al. (1995a) 
A MDS , Driver et al. (1995) 

• HDF , Abraham et al. ( 1 996) 




(c) Late-type Spirals/Irregulars 

-k MDS, Glazebrook etal. (1995a) 
A MDS, Driver etal. (1995) 
• HDF, Abraham etal. (1996) 



(d) All types 




* MDS, Glazebrook etal. (1995a) 
A MDS .Driver etal. (1995) 

• HDF , Abraham et al. (1996) 
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Fig. 5. — Differential number counts for galaxies as a function of apparent magnitude in Isu band. Panel (a), (b), (c) and 
(d) are for ellipticals, early-type spirals, late-type spiral/irregular/vB galaxies and the overall population, respectively. The model 
predictions are shown by lines. 
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Fig. 6. — The evolution of the characteristic number density <j)* of E/SO galaxies with respect to z. <f>* is normalized to unity 
at the present day. 
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Fig. 7. — The evolution of the characteristic number density cj>* of Sdm/Irr and vB galaxies with respect to z, constrained by 
the conservation of the comoving mass density, Eq. [4] . <f>* is normalized to unity at the present day. 




Fig. 8. — Differential number counts for galaxies as a function of apparent magnitude in Hi.e band. Panel (a), (b), (c) and (d) 
are for ellipticals, early-type spirals, late-type spiral/irregular /vB galaxies and the overall population, respectively. The data arc 
taken from NICMOS (Teplitz et al. 1998); the model predictions are shown by lines. 
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(a) 



Metcalfe etal.(1995b) 
Driver etal.(1994) 
Metcalfe et al.(1991) 
Jones etal.(1991) 
Lilly etal.(1991) 
Maddox et al.( 1990b) 
Shanks (1989) 
Tyson (1988) 
Stevenson etal. (1986) 
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a Soiferetal.(1994) 

O HDS Gardner etal.(1993) 

* HMDS Gardner et al.( 1 993) 

• HMWS Gardner etal.(1993) 
□ HWS Gardner etal.(1993) 

* Mobasheretal.(1986) 
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Fig. 9. — Differential number counts as a function of apparent magnitude in B- 
data are exhibited in the figure. Model predictions are shown by lines. 



and K- bands. The sources of observational 
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Fig.lO-(d) 



CFRS data 
17.5 < I AB < 22.5 
completeness=85% 




All 

E/SO 

Sab 

Sbc 

Scd 

Sdm 



A,vB 



Fig.lO-(e) 



CFRS data 
17.5 < I AB < 22.5 
completeness=85% 




B , All 
B , E/SO 
B ,Sab 
B , Sbc 
B , Scd 
B ,Sdm 
B ,vB 



Fig.lO.-(f) 



CFRS data 

17.5 < I AB < 22.5 
completeness=85% 




C, All 

C , E/SO 

C.Sab 

C.Sbc 

C.Scd 

C.Sdm 
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Fig. 10. — The z distributions for galaxies. Fig.lO-(a), (b), and (c) are limited in the magnitude range of 22.5 < B < 24.0, and 
for Scenario A, B, and C, respectively. Fig.lO-(d), (e), and (f) are corresponding to (a), (b), and (c), respectively, but are limited 
in 17.5 < Iab < 22.5. The sources of the observed data are indicated in the figure, and they are shown by histograms. Model 
predictions are indicated by lines, which have been normalized to give the total number of both 2-identified and z-unidentified 
objects. 



